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Abstract 

In the framework of a QCD relativistic potential model we evaluate the form factors de- 
scribing the exclusive decays B — > ntv and B — > K£ + £~. The present calculation extends a 
previous analysis of B meson decays into light vector mesons. We find results in agreement with 
the data, when available, and with the theoretical constraints imposed by the Callan-Treiman 
relation and the infinite heavy quark mass limit. 



The study of the decays 

B -» n£u t (1) 

B -> (2) 

represents a significant part of the experimental programmes at the next proton-proton accelerators 
and at the future 5-factories at SLAC and KEK. The importance of these processes arises from 
the following reasons. The decay (^) allows to measure the product of the Kobayashi-Maskawa 
(KM) matrix element V u b and the form factor describing the decay process^ ; similarly, the decay 
(0) will give access, in appropriate regions of phase space, to the KM matrix element Vj s ; therefore 
these processes would allow to measure fundamental parameters of the Standard Model (SM) of 
the fundamental interactions, to say nothing of the possibility to explore, in both cases, new effects 
beyond the SM. 

It is fair to say, however, that, in spite of the fundamental relevance of the processes (|) and (|D, 
the basic theory of the hadronic interactions, Quantum-Chromo-Dynamics (QCD), is still unable to 
produce clear predictions for the hadronic matrix elements B — > it, B — » K involved in these decays. 
This is due to the lack of a theoretical tool, as powerful as perturbation theory, able to produce 
predictions for the nonperturbative quantities involved in these processes. The most frequently 
used theoretical methods to deal with these problems are based on approximation schemes such 
as lattice QCD or QCD sum rules. These approaches have however their own limitations. In the 
former method the finite lattice size introduces a cut-off in the small momenta, which precludes the 
possibility to make reliable predictions in the small momentum transfer region (Q 2 < 15 GeV 2 ) (for 
recent reviews of lattice QCD predictions for B into light meson transitions see e.g. Q). In the case 
of QCD sum rules or their variant, light cone sum rules, the theoretical uncertainties are dominated 
by the peculiar theoretical tools employed by this method (criteria for stability, hierarchic role of 
the different nonperturbative contributions parametrized by the various condensates) and cannot 
be reduced by adding new terms in the Operator Product Expansion (for a discussion see ||]). 

On the basis of these considerations, in ||] we have presented an analysis of semileptonic and rare 
transitions between the meson B and a light vector meson in a QCD relativistic potential model. 
In [|| we argued that, because of its simplicity, this model might be used as a viable alternative to 
the more fundamental, but still limited theoretical approaches we have discussed above. It is the 
aim of this paper to extend this analysis to the decays (H) and @. 

To begin with, we review the main features of the QCD relativistic potential model. It is 
a potential model because the mesons are described as bound states of constituent quarks and 
antiquarks tied by an instantaneous potential V(r). It is a QCD model because the potential is 
modelled according to the theory of the hadronic interactions, i.e. it has a confining linear behaviour 
at large interquark distances r and a Coulombic behaviour ~ —a s (r)/r at small distances, with 
a s (r) the running strong coupling constant: in practice the interpolating Richardson's potential 
V(r) is used 0, cut-off at very small distances (of the order of the inverse heavy meson mass) 

"There is one form factor contributing to ([!]) for a massless lepton. 
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to take care of unphysical singularities introduced by the relativistic kinematics ||. Finally it is 
a relativistic model because the wave equation used to obtain the meson wave function ^ is the 
Salpeter equation embodying the relativistic kinematics: 



V 2 + mf + J -V 2 + m 2 , + V(r 



#(r) = Ml'(r) , (3) 



where, for heavy mesons made up by a heavy quark Q and a light antiquark, 1 refers the heavy 
quark and 2 to the light antiquark. The relativistic kinematics plays an important role when 
at least one of the two quarks constituting the meson is light, as in our case, and represents an 
improvement in comparison with the approach based on the non-relativistic quark model. In (|3|) 
M is the heavy meson mass that is obtained by fitting the various parameters of the model, in 
particular the b-quark mass, that is fitted to the value mj, = 4890 MeV, and the light quark masses 
m u ~ = 38 MeV, m s = 115 MeV f\ . The S-meson wave function in its rest frame is obtained 
by solving (|3|); a useful representation in the momentum space was obtained in || and is as follows 



ip(k) = Wmsa 3 e~ ak , (4) 

with a = 2.4 GeV -1 and k = \k\ the quark momentum in the B rest frame. 

The constituent quark picture used in the model is well suited for the mesons comprising at least 
a heavy quark; for light mesons other dynamical features, not accounted for by this simple picture, 
should be incorporated, e.g. the nature of pseudo Nambu-Goldstone bosons of 7r's and -fT's and the 
presence of important spin — spin terms in V(r), not included in the Richardson's potential (their 
neglect for heavy mesons is justified by the spin symmetry of the Heavy Quark Effective Theory 
(HQET)Q which is valid in the limit niQ — > oo). The solution adopted in H was to avoid, for 
light mesons, the constituent quark picture and to describe their couplings to the quark degrees of 
freedom by effective vertices. This assumption produces a set of rules that are used to compute 
the quark loop of fig. [lj, i.e. the diagram by which the hadronic amplitudes describing the decays 
d]]) and (H) are evaluated. They are as follows. 

1) For a light pseudoscalar meson M (= tt , K) of momentum p' we write the coupling 

^'-75, (5) 



M 



where Jm = fir = 130 MeV or Jm = Jk = 160 MeV . The normalization factors N q , N q i for the 
quark coupled to the meson are discussed below. 

2) For the heavy meson B in the initial state one introduces the matrix: 



b = ±mJ mqmb ^±^ H 75) ^L±^ (6) 

V3 V m q m b + qi ■ q 2 2m b 2m q 



^Data on the heavy meson spectra are not of great help in fitting light quark masses, which, therefore, are not 
accurately determined in the model; its predictions, however, are not sensitive to m u , m,d, m s values in most of the 
available kinematical range. 
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where m b and m q are the heavy and light quark masses, , q% their 4— momenta. The normalization 

f d?k 

factor corresponds to the normalization < B\B >= 2 tub and J — ^ 3 \ip(k)\ 2 = 2itlb already 

embodied in (0). One assumes that the 4— momentum is conserved at the vertex Bqb, i.e. q^+q^ = 
= B meson 4— momentum. Therefore q± = (E b , k), q% = (E q , —k) and 



E b + E q = m B . (7) 

3) To take into account the off-shell effects due to the quarks interacting in the meson, one introduces 
running quark masses m(k), to enforce the condition 



E=\]m 2 {k) + \k\ 2 (8) 

for the constituent quarks. For the kinematics of the decays (|l|) and it is sufficient to introduce 
the running mass only for the heavy quark j] 

m b = m b (k) , (9) 

defined by the condition 



1 ml + \k\ 2 + \Jm 2 b + \k\ 2 = m B . (10) 
4) The condition m 2 > implies the constraint 



2 2 

< k < k M = — | 1 , (11) 

ZrriB 



on the integration over the loop momentum k 

d 3 k 



(12) 



(2vr) 3 ' 

5) For each quark line with momentum q and not representing a constituent quark one introduces 
the factor 

x G(q 2 ) , (13) 

4 - "v 

where G(q 2 ) is a shape function that modifies the free propagation of the quark of mass m q i in the 
hadronic matter. The shape function 

2 2 

„ m r — m„, 

G(Q 2 ) = — 2 1 ( 14 ) 

rriQ — q z 

was adopted in |H] ; the value of the mass parameter vtlq was determined in Q by the experimental 
data on the B — > K*7 decay. A range [1.2, 7.6] GeV 2 of possible values of m 2 G was obtained. 

6) For the hadronic current in fig. |l] one puts the factor 

N q N q ,F» , (15) 

*By this choice, the average < m,b(k) > does not differ significantly from the value m;, fitted from the spectrum, 
see Hi for details. 
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where is a 4 x 4 matrix. We shall consider T M = 7^ and = a^ u q v (with cr^ = i/2 [7^,7"]) 

(if q = constituent quark) 
1 (otherwise) . 



The normalization factor N q is as follows 




(16) 



7) For each quark loop one puts a colour factor of 3 and performs a trace over Dirac matrices. 

This set of rules can now be applied to the evaluation of the matrix element < M(p')\q'T^b\B(p) > 
with the result: 



< M(p')\qT^b\B(p) >=Vsj ^0[k M ~ A#(*0y/ 



m q m b 



m q m b + qi-q 2 



Tr 



2mi 



2m n 



?1 



• (17) 



From this expression one can obtain the relevant formulae for the various form factors. With 
q = p — p', we write 

< M{p')\q'^b\B(p) > 



Fi(q 2 )(p+pT + ^ 2 mj V (F (q 2 ) - 



(18) 



< M(p')\ i q'a^q u b\B(p) > 



Mq 2 



rriB + rriM 



{p + p'Yq 2 - (m 2 B - m 2 M )q» 



where 



Fi(q 2 ) = f + {q 2 ) 
Fo(q 2 ) = U(q 2 ) + 



m%-m 2 M 



f-(q 2 



(19) 



(20) 



In ( |j~8| ) and (19) we shall consider M = ir or M = K since both cases are of physical interest if we 
wish to consider not only semileptonic and radiative transitions, but also nonleptonic decays. 



The calculation of the trace and the integral in (|17| ) is straightforward and is similar to the 
one obtained in || for B —* p, B —* K* transitions. For all the form factors we write F(q 2 ) = 
F(q 2 ,m q /) — F(q 2 jUic), where, for the various form factors, we have 



F (q 2 ,x) 



^ 2 h 



dk k 2 ip(k) 



M (m B - m M ) Jo ^E q E b [m 2 B - {m b - m q ) 2 } 



dz 



-1 m M — 2E q (mB — q°) + rn 2 — x 2 + 2\q\ k z 
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[2E q (m B - q c 



m 



M 



2\q\ k z][(m q E b + m b E q )q° + (m b - m q ) \q\ k z] + 



m a + m 



(q 2 - m B q°) m B - (m b - m q ) 2 + 2 m 2 B \q\ k z 



(21) 



Fi{q 2 ,x) 



\/6 f k M 



dk k 2 if;(k) 



n fa Jo JE q E b [m B -{m b -m q y\ 



dz 



m 2 M - 2E q (niB - q°) + m 2 - x 2 + 2\q\ k z 



[2E q (m B - q°) - m 2 M - 2\q\ k z] 



\q\(m q E b + m b E q ) + q° k z (m b - m q ) 
m B \q\ 



m q + niqi 2 k z {tub q° — q 2 ) — \q\ [m B — {m b — m q ) 2 ] 
2 W\ 



+ 



(22) 



h{q\x) 



V6 tub + rriM f kM 
o 



i 



dk k 2 ^(k) 



E q E b [m B - (m b - m q ) 2 ] 



dz 



-l m 2 M — 2E q (niB — q°) + m 2 q — x 2 + 2\q\ k z 
[2E q {m B - q°) - m 2 M - 2\q\ k z] ^ + 



(m q + m ql ) 



\q\ 

-(m b E q + m q E b ) \q\ + (m B - q°) k z (m b - m q 



m B \q\ 

In these equations qo is the time component of four- momentum q^, 

z = cos(9) , 



(23) 



(24) 



with 9 the angle between k and the direction of transferred momentum q. We note that, for M = tt, 
Trig = m q i = m u and fu = fir, while, for M = K, m q = m u , m q i = m s and Jm = fx- 

Before discussing our numerical results in detail let us compute i ? o(9 2 ) f° r q 2 = m e — m Mi ^ n 
the chiral limit Fo(m B — m 2 M ) ~ Fo(m B ) must obey the Callan — Treiman relation 



7-i ( 2 \ fB 

F o(m B ) = -j- 

Jit 



(25) 



This is therefore a consistency test to be satisfied by the model. We have numerically evaluated 
FQ 7T {m 2 B ) for different values of the parameter vug an d we have obtained the result F^ n (m B ) ~ 
1.48, almost independent of mo- This result should be compared to f B / fn — 1.58, which is obtained 
using f B = 0.2 GeV, i.e. the value computed in || using the present model. The small discrepancy 
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in the Callan-Treiman relation (~ 6%) may be attributed to the deviations induced in the B meson 
wave function by the chiral limit that are not accounted for by this calculation. We expect however 
that these differences vanish if, in addition to the chiral limit, one also takes the infinite heavy 
quark mass limit; as a matter of fact one can verify rather easily, using the previous formula for Fq 
and the expression in || for fs, that the Callan-Treiman exactly holds in the combined mj — > oo 
and rriM —> limit. 

Let us now consider the form factor F\{q 2 ) (respectively /t(<? 2 ))- Our numerical results for 
the central value of rog, i.e. mg = 1.77 GeV, show that the q 2 — behaviour of this form factor 
is increasing (resp. decreasing) for both small and moderate values of q 2 , independently of the 
value of the mass parameter mg introduced in eq. (|l3|) . This behaviour should hold also at large 
q 2 (g 2 > 15 GeV 2 ) due to the effect, in this region, of a pole in the q 2 functional dependence, 
predicted by the dispersion relation. Differently from our analysis of Fo(q 2 ), we cannot pretend, 
however, to extend the validity of our predictions for F\{q 2 ) and Jt(q 2 ) at the extreme values of 
q . The difference between the two cases is as follows. In the case of F , the pole (with J p = 0+) 
contribution to this form factor vanishes in the chiral limit and has therefore a minor impact on the 
q 2 behaviour. On the contrary the form factors F\ and /t, have a non vanishing polar contribution 
which becomes larger and larger with increasing q 2 . While we expect that this behaviour become 
visible well before the pole, at extreme values of q 2 the diverging behaviour induced by such a 
contribution cannot be reproduced by the model. As a matter of fact, for larger values of q 2 , \q[ 
becomes smaller and smaller, and, therefore, the model becomes sensitive to the actual values of 
the parameters, in particular the light quark masses that are not accurately fitted by the available 
experimental data (see above). Therefore we can consider that our predictions are reliable in the 
range (0, 15) GeV 2 ; at q 2 = we get 

Ff n (0) = Ff n (0) = 0.37 ± 0.12 

F^ K (0) = F^ K (0) = 0.26 ± 0.08 

/^(0) = -0.14 ±0.02 

f BK {Q) = -0 m +O-OB (26) 

The central values are obtained for mg = 1.77 GeV, which is the best fit of the parameter mg 
found in || by the experimental branching ratio B(B — ► K #r y), whereas the theoretical uncertainty 
is obtained by varying mg in the range [1.1, 2.8] GeV. The results for B — > tt refer to charged 
pions. 

Let us now consider the g 2 -behaviour of the form factors. We introduce the two-parameter 
function for the three form factors 

nf) = ; (27) 

1 - 

here ap ,bp are parameters to be fitted by means of the numerical analysis and F(0) is given in 
eq. (EB); to allow a comparison with other approaches we perform the analysis up to q 2 = 15 GeV 2 , 
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both for 7r and K mesons. We collect the fitted values in table [l] and report the g 2 -dependence in 

fig- B I- 





F(0) 


ap 


b F 


F(0) 


ap 


bp 




F B* 


0.37 


0.60 


0.065 


0.26 


0.50 


0.39 


pBK 




0.37 


1.1 


0.44 


0.26 


1.2 


0.56 


zpBK 


fB-K 

It 


-0.14 


0.92 


0.21 


-0.09 


0.76 


0.76 


f BK 

It 



Table 1: Parameters appearing in eq. (27) for different B form factors. 



From table |T| and from fig. |], |3|, one can see that Fo(q 2 ), F^ n (q 2 ) and fp n {q 2 ) have a q 2 behaviour 
similar to a single pole. For F^ K (q 2 ) and fp K (q 2 ) there are significant deviations from this 
behaviour. We do not have yet experimental data to test these predictions and we shall limit to 
compare our results with other theoretical approaches; before doing that, let us discuss the infinite 
heavy quark mass limit of the model. In the framework of the Heavy Quark Effective Theory, which 
corresponds to the m& — * oo limit, there is a constraint to be satisfied by the three form factors, 
i.e., the relation originally found in 0: 



fr{q ) = 



F(a 2 \ (m 2 m 2 \ F 0(Q 2 )-Fl(q 2 ) 
l<i(q )- [m B - m M — 



q 2 



(28) 



2 TflTj 

This relation holds in the limit m;, — > oo and for high q 2 (q 2 ~ q^nax)- We have checked that this 
relation formally holds in our model in the infinite heavy quark mass limit and for q 2 ~ q^ax- F° r 
the actual value of mp (= 5.28 GeV) and for the transition B — > it the situation is as follows. We 
choose q 2 M = 15 GeV 2 , i.e. the maximum value at which we can trust our predictions and we find 
numerically F 1 (q 2 M ) ~ 0.54, F (qj d ) ~ 0.71, fT{q 2 M ) ^ -0.27. Therefore the relation (||) has a 
significant violation of 50%, that may be attributed to the fact that we are still far from g^ ai , and 
0(l/mb) corrections are large. Similar results are obtained for the B — > K transition. 

Let us finally comment on the scaling laws of the form factors at large q 2 that can be helpful in 
using the heavy flavour symmetry to relate the form factors of B and D mesons [ 10 1 . From eqs.(^2|) 
and (|2^) the following behaviours can be formally derived: 

Fi{q 2 max ) ~ V^b (29) 

hiqlnax) ~ \/m> ; (30) 

they are in agreement with the pole (vector meson) dominance of the form factors observed in || 
[ |10| ; moreover in the chiral limit, from eq.(^Tj) one gets: 

F (q 2 max ) « -L . (31) 



Let us now compare our work with other theoretical approaches. In table |2| we compare our 
outcome for the values at q 2 = with the results of QCD sum rules and lattice QCD calculations 
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(for other work on this subject see, e.g. We observe that our results are in agreement, within 

the theoretical uncertainties, with the determinations obtained by light cone sum rules (LCSR) 
0, lattice |] and lattice + LCSR 





This work 


LCSR 1 12] 


SR [13| 


Latt. [§ 


Latt. + 
LCSR [14| 


Ff-(O) 
f^(0) 


0.37 ±0.12 
-0.14 ±0.02 


0.30 ± 0.04 
-0.19 ±0.02 


0.24 


0.27 ±0.11 


0.27 ±0.11 


Ff K (0) 
f* K (0) 


0.26 ±0.08 
-0 0Q+ - 05 

u - uy -0.02 


0.35 ±0.05 
-0.15 ±0.02 


0.25 
-0.14 







Table 2: Comparison of the results coming from different approach to evaluate form factors. 



As for the q 2 dependence, we have not reported the predictions of other theoretical approaches, 
since they qualitatively agree with our calculations. In absence of detailed experimental data on the 
form factors, the best we can do to test the model is to use data on the partial width T(B — > irlv). 
To perform this comparison we must, however, extrapolate the c/ 2 -behaviour obtained by eq. ( p7[ ) 
and tables || and [2], and valid in the region (0, 15) GeV 2 , to the whole q 2 range. This procedure 
implies an uncertainty which is difficult to assess, but should not be extremely large due to the 
phase space limitation at high q 2 . We obtain 

BTZ(B° tt+Tp) = 1.03 ( |T4fe| 10~ 4 , (32) 

\3.2 10 6 J 

to be compared to the experimental value B1Z(B — ► n£u) exp _ = (1.8 ± 0.6) 10~ 4 JT5[ ] . Therefore our 
result is compatible with the present range of the KM matrix element V u b = [1.8, 4.5] x 10~ 3 ; 
the preferred range of values selected by the model and by the present experimental limits on V u b 
is V ub = (4.0 ±0.5) x 10~ 3 . 

We conclude our analysis by summarizing our results. We have used a QCD relativistic potential 
model, introduced in Q, to study the weak and radiative transitions B — > n, K. We have computed 
the relevant form factors and tested the Callan-Treiman and the Isgur-Wise relation. The former 
relation, valid in the chiral limit, is satisfied at the 6% level, while the latter, valid in the nib 00 
limit, has significant violations, due to 0(1/ nib) corrections. Our result for the branching ratio 
BTZ(B — » ttIv) agrees with the experimental data. 
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stage of this work. 
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-Q2 

Figure 1: Quark loop diagram describing the matrix element < M(p')\J fl \B(p) >; M is a light 
pseudoscalar meson, = q'T^b is the current inducing the decay and is a combination of Dirac 
matrices. 
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